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Abstract. We have used the SEST 15-metre and Onsala 20-metre telescopes to perform deep (r.m.s. ^ 30 mjy) 
integrations of various molecular rotational transitions towards damped Lyman-alpha absorption systems (DLAs) 
known to occult millimetre-loud quasars. We have observed 6 new systems and improved the existing limits for 
11 transitions. These limits may be approaching the sensitivities required to detect new systems and we present 
a small number of candidate systems which we believe warrant further observation. 
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1. Introduction 

Millimetre-band molecular absorption systems along the 
line-of-sight toward quasars can provide a powerful probe 
of cold gas in the early Universe. Wiklind & Combes (1994; 
1995; 1996b) have used such absorption lines to study a 
variety of properties of the absorbers themselves (e.g. rela- 
tive column densities, kinetic and excitation temperatures, 
filling factors). Constraints on the cosmic microwave back- 
ground temperature can also be obtained by comparing 
the optical depths of different rotational transitions (e.g. 
Wiklind & Combes 1996a). If the background quasar is 
gravitationally lensed, time delay studies can yield con- 
straints on the Hubble constant (e.g. Wiklind & Combes 
2001). However, these studies have so far been limited 
by the paucity of mm-band molecular absorbers: only 4 
such systems are currently known-towards TXS 0218-f 357 
(Wiklind & Combes 1995), PKS 1413-fl35 (Wiklind & 
Combes 1997), TXS 1504-^377 (Wikhnd & Combes 1996c) 
and PKS 1830-211 (Wikhnd & Combes 1998). 

Recent attention has focused on using molecular ab- 
sorption lines as a probe for possible changes in the fun- 
damental constants on cosmological time-scales. Detailed 
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studies of the relative positions of heavy element opti- 
cal transitions compared with laboratory spectra favour 
a smaller fine structure constant (a = e'^ /he) at redshifts 
0.5 < z < 3.5 at the 4.1(t significance level (Murphy et al. 
2001b, Webb et al. 2001). The observed fractional change 
in a [Aa/a — (—0.72 ±0.18) x 10"^] is very small, so sys- 
tematic errors have to be carefully considered. However, a 
thorough search for systematics has not revealed a simpler 
explanation of the optical results (Murphy et al. 2001c) 
and so independent constraints at similar redshifts are re- 
quired. 

Comparison of molecular rotational (i.e. mm-band) 
and corresponding Hi21-cm absorption line frequencies 
has the potential to constrain changes in a with a frac- 
tional precision ~ 10~^ per absorption system. The ratio 
of the hyperfine (H i) transition frequency to that of a 
molecular rotational line is oc y = a^gp and so any varia- 
tion in this will be observed as a difference in the apparent 
redshifts of these lines (Drinkwater et al. 1998). 

From spectra of two of the known high redshift ab- 
sorbers, PKS 1413-1-135 and TXS 0218+357, CariUi et al. 
(2000) and Murphy et al. (2001a) have obtained con- 
straints on Ay/y consistent with zero variation at red- 
shifts Zabs — 0.6847 and 0.24671, respectively. However, 
the major uncertainty in the mm/Hi comparison is that 
intrinsic velocity differences between the mm and H i ab- 
sorption lines can be introduced if the lines-of-sight to the 
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Table 1. Top: The northern sample. The radio flux densities at various GHz frequencies of the background quasars are 
given (see Curran et al. 2002 for details). Note that a "w" denotes a variable flux density. 2;abs is the DLA redshift and 
the final column gives the molecules most commonly detected in the four known absorbers (see Wiklind & Combes 
references), which fall into the Onsala 84-116 GHz band. (*HCO+ ^ 1 in B 08279+5255 and *CS 2 ^ 3 in B 
1017+1055 are 7 mm observations). Note that CS has so far only been detected in PKS 1830-211 (Wiklind & Combes 
1996c) but since none of the commonly detected molecules fell into the band at the DLA redshift we timed to this 
dense gas tracer. Bottom: The southern sample. The final column lists the molecules most commonly detected which 
fall into the SEST 78-116 GHz and 128-170 GHz bands. 
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Table 2. Molecular lines previously searched for in radio-illuminated DLAs. All of the observations were performed 
with the IRAM 30-m telescope, except for HCO+ 2^3 towards B 0458-020 which was observed with the SEST. 
Again, apart from those marked "*" which are from Wiklind & Combes (1995), the radio flux densities are from the 
various sources cited in Curran et al. (2002). cr is the lowest r.m.s. noise obtained for this frequency (with the reference 
given) recalculated for a channel width of 10 km s~^(see footnote 3). Note that Wikhnd & Combes (1994b) is a search 
for CO emission, although any absorption features at the appropriate redshift should be apparent. 
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millimetre wave and radio continuum emission regions of 
the quasar differ, as is certainly the case for PKS 1413+135 
and TXS 0218+357 (CariUi et al. 2000). Thus, a statisti- 
cal sample of mm/Hi comparisons is required to indepen- 
dently check the optical results. 

One systematic approach to finding new high redshift 
molecular absorbers is to target high column density sys- 
tems with known redshifts. A convenient sample is there- 
fore the damped Lyman-alpha absorbers defined to have 
neutral hydrogen column densities A'hi ^ 10^" cm~^. 
We therefore compiled a catalogue of all known DLAs 
and shortlisted those which arc illuminated by radio-loud 
quasars (Curran et al. 2002)[| Of these, seven in the north- 
ern sky have measured millimetre fluxes, and in the south 
there are around a dozen DLAs illuminated by millimetre- 
loud quasars, of which we observed the four loudest (Table 
0). 

We note that Wiklind & Combes (1994b; 1995; 1996b) 
searched with null results for redshifted molecular emis- 
sion and absorption towards 12 and 46 quasars, respec- 
tively. However, only 11 are occulted by DLAs and not all 
of these are radio- loud (Table ^ . This motivates a more 
systematic search for millimetre absorption systems asso- 
ciated with DLAs. In this paper we present the results of 
our first search: the DLAs which occult known millimetre- 
loud quasars with the SEST^ and Onsala 20-m telescopes. 

2. Observations 

The northern 3 mm and 7 mm observations were per- 
formed in February and April 2002, respectively, with 
the 20-m telescope at Onsala Space Observatory, Sweden. 
For the 3 mm observations the SIS receiver was tuned 
to single-sideband mode and the backend was a hybrid 
correlator with a bandwidth of 1280 MHz and a channel 
separation of 0.8 MHz. We used dual-beam switching with 
a throw of 12' in azimuth, and pointing errors were typi- 
cally 3" r.m.s. on each axis. The intensity was calibrated 
using the chopper-wheel method and typical system tem- 
peratures, on the T^-scale, were around 300 K. For the 
7 mm observations we used the 43 GHz SIS receiver and 
again the backend was the hybrid correlator with a band- 
width of 640 MHz and a channel separation of 0.4 MHz. 
Since there is no beam switch capability at these frequen- 
cies, we removed the ripple caused by standing waves by 
subtracting the Fourier components in a transform of the 
spectra. Typical system temperatures, on the T^-scale, 
were around 200 K. 

The southern sources in the sample were observed in 
April 2002 with the 15-m SEST at La Silla, Chile, using 
the 100 GHz and 150 GHz SESIS receivers. These were 
tuned to single-sideband mode and typical system temper- 
atures, on the T^-scale, were ~ 200 K for the RXlOO and 
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^ Available from 

^ The Swedish-ESO Sub-millimetre Telescope is operated 
jointly by ESO and the Swedish National Facility for Radio 
Astronomy, Onsala Space Observatory, Chalmers University 
of Technology. 



« 250 K for the RX150. The backends were acousto-optic 
spectrometers with 1440 channels and a channel width 
of 0.7 MHz. As with the Onsala 3 mm observations, we 
used dual-beam switching with a throw of about 12' in 
azimuth, and pointing errors were typically 2" r.m.s. on 
each axis. Again, the intensity was calibrated using the 
chopper-wheel method. 

3. Results and Discussion 

Upon the removal of a low order baseline and subsequent 
averaging of the data for each quasar, no absorption fea- 
tures of > 3cr/channel were foundf} In Table | we sum- 
marise the derived upper limits together with the previ- 
ously published results. Comparing the optical depth lim- 
its with those in the four known absorbers, the strongest 
absorber, towards PKS 1830-211, has Tmm ~ 1 for HCO+ 
1 ^ 2 and Tnun ~ 0.5 for HCO+ 2^3. The weakest, to- 
wards PKS 1413+135, has r^m ~ 0.1 for both CO ^ 1 
and HCO^ 2^3. The two remaining systems have op- 
tical depths similar to that towards PKS 1830-211 and 
so our new observations and the previous surveys should 
have been sensitive to (at least) CO and HCO"*" at similar 
strengths to those in 3 of the 4 known systems. Only for 
the observations of B 0235+1624, HCO+ 2^3 towards 
B 2136+141 (IRAM) and CO ^ 1 towards B 1451-375 
are the surveys sensitive enough, particularly for resolu- 
tions of ^ 10 km s~^, to detect an absorption system of 
a strength similar to that towards PKS 1413+135 (See 
Table 4). 

From the optical depths we may estimate column den- 
sity limits for each transition. In Table ^ we give the 3a 
optical depth limits according to a velocity resolution of 
10 km s~^ as well as at a resolution of 1 km s~^. We give 
the former in order to show the limits corresponding to a 
visual inspection of the spectra^ and facilitate a more di- 
rect comparison with the previous surveys. By multiplying 
the latter with the expected width of a line, the velocity 
integrated optical depth for a 3cr detection at At; = 1 
km is obtained. For all of the limits, assuming LTE 
conditions, we can estimate the total column density of 
each transition from 



rdv, 



(1) 



where is the rest frequency of the J — > J + 1 transi- 
tion, and Afj, are the statistical weight and the Einstein 
A-coefficient of the transition, J rdv is the velocity inte- 
grated optical depth and / is the product of the partition 
function for an excitation temperature, (assumed w 10 
K), and e^'/'='^-/(l - e-^"/''^-) (see Wiklind & Combes 
1995, 1996b, 1999 for details). From the optical depth lim- 
its, we derive limits on the total column density per unit 
velocity, thus normalising the limits which depend on the 
spectral resolution to which the data have been smoothed. 



^ See |ittp://www.phys.unsw.edu.au/'-^sjc/dla-figl.ps.gz for 
the spectra and corresponding r.m.s. noise levels. 
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Table 3. Summary of our (Table |l]) and the previously published results (Table |]). F is the visual magnitude of 
the background quasar, A^hi [cm^^] is the DLA column density from the Lyman-alpha line and T21 cm is the optical 
depth of the redshifted 21 cm Hi line (see Curran et al. 2002). The optical depth of the relevant millimetre line is 
calculated from r = — ln(l — 3(Trnis/'S'cont)i where CTi-ms is the r.m.s. noise level at a given resolution and Scant is the 
continuum flux density, estimated from the values at the neighbouring frequencies (Tables]^ and|^). This is done for 
resolutions of 10 km s^^ [T,„m(10)] and 1 km s~^ (rmm), where for the latter we quote only the best existing limit. 
For all optical depths, 3a upper limits are quoted and "-" designates where 3a > Scoot, thus not giving a meaningful 
value for this limit. Blanks in the T21 cm held signify that there are no published Hi absorption data for these DLAs. 
The final column gives the best existing limit of the column density per unit velocity estimated for the transition 
[cm~^(km s~^)~^] (see main text). 
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Notes: Where we have observed towards the same quasar *denotes Wiklind & Combes results (see Table g for details). "H I ab- 
sorption a.t z = 0.52385 (Briggs & Wolfe 1983). ''Hi absorption at 2 = 2.03945 (Wolfe et al. 1985; Briggs et al. 1989). Assumed 
flat spectrum, i.e. 545 = 5*87 — S350 ~ 0.08 Jy. Assumed 5*44 « 0.015 Jy and 5*102 ~ 0.008 Jy; "^Sus ~ 0.011 Jy, based on the 
two available measured values (Table 111) . 



Since we have only upper limits, and thus no knowl- 
edge of the width of any line which may be hidden in the 
noise, as well as uncertainties in the conversion to H2 col- 
umn densities for molecules other than CO, it is difficult 
to draw meaningful comparisons between iVfji and iVmm- 
However, by assuming an absorption line of FWHM 20 
km s~^ (as in the case of 3 of the 4 known absorbers), we 
estimate for the lowest optical depth limit (CO 1 — > 2 in B 
0235+1624) a value of iVco ~ 6x 10^''^ cm'^ at the 3a level 



and a resolution of 1 km s"^, i.e. Nco/Nm ^ 2 x 10"^. 
According to A'^h2 lO^'A^co (e.g. Wiklind & Combes 
1998), the molecular to atomic hydrogen column density 
ratio is A^Ha ~ 2% A^hi, which is similar to the lower 
limit estimated from the optical detection of molecular 
hydrogen^ in the DLA towards B 0528-2505 (CariUi et 

^ In the case of z > 1.8 sources, the ultra-violet lines of H2 are 
redshifted into the optical window, making molecular hydrogen 
readily observable at these frequencies. As well as towards B 
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al. 1996). That is, although we have no clear detections 
of any molecular absorption lines, the current limits may 
be approaching those sufficiently low in order to detect 
molecular absorption in damped Lyman-alpha systems. 

Bearing this in mind, we then analysed the spectra 
for tentative features within a range corresponding to the 
uncertainty in the DLA redshift^. In each case the con- 
tinuum level was defined using a polynomial baseline fit 
and for each channel we generated a Icr error from the 
r.m.s noise in a window of 2N^in + 1 channels centered on 
that channel. Using this error array, we identified absorp- 
tion features as series of n channels over which a deviation 
from the continuum level was observed with significance 
> (Tiiin standard deviations. Note that the significance of 
an absorption feature, cr„, is an overestimate because the 
flux in adjacent channels is (positively) correlated [each 
spectral resolution element is sampled at approximately 
the Nyquist rate, see Murphy, Curran & Webb (2002) for 
details] . 

According to the "expected" (5los = 100 km s^^ veloc- 
ity differences in the optical and millimetre wave lines-of- 
sight^, the only candidate features identified in our anal- 
ysis are HCO+ 2 ^ 3 and CO 3 ^ 4 towards B 0458-020 
as well as the HCO+ 2^3 line towards B 0834-201. 
As seen in Table 4, however, increasing the line-of-sight 
velocity difference to (Slos — 200 km s~^ significantly 
increases the number of candidates. This suggests that 
either the authenticity of the candidate features is ques- 
tionable and/or we have underestimated the uncertainties 
in the DLA redshifts. It is interesting to note that all but 
one of the candidate lines occur in DLAs with high col- 
umn densities (A^hi <^ lO'^^ cm^'^). The one exception is B 
0834-201 which has the highest visual magnitude of these 
candidates. We emphasise that Table 4 is only intended as 



0528-2505 (Foltz, Chaffee & Black 1988) molecular hydrogen 
has also been detected in the DLAs occulting the radio-quiet 
quasars B 0000-2620 (Levshakov et al. 2000), B 0013-0029 
(Ge & Bechtold 1997; Petitjean, Srianand & Ledoux 2002), 
B 0347-3819 (Levshakov et al. 2002), B 1232+0815 (Ge & 
Bechtold 1997; Srianand, Petitjean & Ledoux 2000) and the in- 
ferred (Wolfe et al. 1995) DLA towards B 0551-3637 (Ledoux, 
Srianand & Petitjean 2002). 

^ Although we have searched for possible absorption lines 
close to the expected redshifts, it may possible that the line- 
of-sight to the mm-band and optical emission regions of the 
quasars may differ ((5los > 0), giving rise to intrinsic off- 
sets between any millimetre absorption and the known DLA 
redshifts. Drinkwater et al. (1998) compared millimetre and 
Hl21-cm Galactic absorption profiles and found that the in- 
dividual velocity components in the profiles corresponded to 
within (5los = 1.2 km s"'^. Carilli et al. (2000) argued that 
Slos ~ 10 km s~^ is typical of the velocity dispersion of the 
interstellar medium in galaxies, and that it may even be as 
high as Shos ~ 100 km s~^. In order to account for this, 
we searched for possible absorption lines in a velocity region 
[— (5los — cTDLAjf^LOS + ctdla] Centered on the expected DLA 
frequency with (5los ~ 100 km s~^ and ctdla the uncertainty 
in the DLA redshift. 



a shortlist for follow up observations, which may or may 
not yield absorption systems at the listed redshifts. 

Of all the candidate systems, the most promising 
and visually striking example of a possible absorption 
line is the 151.71 GHz feature towards B 0458-020 (Fig. 
|l|), which has a high visual magnitude as well as the 
highest neutral hydrogen column density of the sample. 
Comparing this "detection" of -8 mK (0.24 Jy at SEST) 



0458-020 




1 1 11 III 1 


C03-4 


DLA HI 

- 1 , 1 i 1 









I 151.4 151.6 151.8 152 



Frequency [GHz] 

Fig. 1. CO 3 ^ 4 at z = 2.039 towards B 0458-020 at 
the full spectral resolution of 1.38 km s^^. The antenna 
temperature of the peak channel is 8 mK (c.f. cr = 3 
mK) and the integrated main-beam intensity of the line 
is —0.10 ± 0.02 K km s~^. The expected frequencies of 
the line for the (approximate) DLA redshift of 2.0399 and 
Hi redshift of 2.03945 (Turnshek et al. 1989, Wolfe et al. 
1995) are shown. 

with the estimated flux density at 152 GHz (0.62 Jy, esti- 
mated from the neighbouring values in Table 0) gives an 
optical depth of Tmm « 0.5. The r:; 9 km s"^ FWHM of 
the line gives a column density estimate of A^co ~ 8 x 10^^ 
cm-2 or A^Ha ~ 10^^ cm'^, i.e. ~ 20% of A^hi- This com- 
pares well with the column densities and their ratios for 
the known absorbers, i.e. 10 — 40% (Wiklind & Combes 
1994, 1995, 1996c, 1999, CariUi et al. 1998), although it is 
considerably higher than the rest of the DLA sample. 

4. Summary 

We have observed 18 transitions for molecular absorption 
over the redshift range 0.2 > z <; 3 in 11 damped Lyman- 
alpha absorption systems which lie along the line-of-sight 
towards radio-loud quasars. Of these, 9 quasars have flux 
densities ^ 0.1 Jy in the millimetre band, making our 
observations sensitive to optical depths ^ 0.1 at a 3(t level 
and a spectral resolution of 10 km s~^, thus improving the 
limits for 11 transitions over the previous results. 

From these observations and the previously published 
results, we estimate an upper limit of sa 2% for the molec- 
ular hydrogen to neutral atomic hydrogen column density 
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Table 4. Candidate molecular absorption systems located within « ±200 km s^^ of the DLA redshift. Vobs is the 
observed frequency, cr„ is the significance of the putative feature, Zmm is the absorption redshift and Z21 cm is the 
redshift of the Hi line (see Curran et al. 2002). The final column gives the offset in km between the z^ni and the 
nominal value of zdlAj with our estimate of uncertainty in the DLA redshift"^ quoted. 



DLA 


Transition 


Uobs [GHz] 






2:21 cm 




B 0248+430 


CS 2 3 


105.3701 


4.2 


0.39479 


0.3941 


190 ± 20 


B 0458-020 


HCO+ 2^3 


87.9994 


3.8 


2.04045 


2.03945 


50± 10 






88.0229 


3.2 


2.03964 




-30 ± 10 




CO 3 ^ 4 


151.7144 


4.3 


2.03888 




-100 + 10 


B 0738+313 


CO ^ 1 


94.3116 


4.4 


0.22224 


0.2212 


240 + 20 


B 0834-201 


HCO+ 2^3 


98.5135 


4.4 


1.71595 




100 + 100 


PKS 1413+135 


HCO+ 2^3 


214.6110 


4.6 


0.24671 


0.24671 





Note: The last entry is a known absorption system included as a guide. A spectrum was obtained after ~ 30 hours with SEST 
(August 2001) which may well have been rejected upon visual inspection^. 



ratio. This is the lower limit derived by Carilli et al. (1996) 
based on a molecular hydrogen detection and a redshifted 
Hi 21 cm non-detection. This is consistent with the re- 
sults of Liszt (2002) who suggests that the low metallicity 
of the Lyman-alpha system is not favourable for the for- 
mation of molecules, but our limit perhaps suggests that 
we are close to the values necessary to detect absorption 
lines in at least the high column density (A'hi ^ 4 x 10^-'^ 
cm~^) DLAs. Subsequently, from an analysis of the spec- 
tra we find several significant features located close to the 
expected redshifts, which do, however, require confirma- 
tion. Therefore follow up observations of the candidates 
given in this article should be the next step in searching 
for the elusive high redshift molecular absorption systems. 
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